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Recent experiments and theories suggest that spin-orbit couphng (SOC) can 
lead to new phases of quantum matter with highly non-trivial collective quantum 
effects. Two such phases are the quantum spin Hall insulator and the 



strong topological insulator (STI) 3, ^, 3] both realized in the vicinity of a 
Dirac point {s] but yet quite distinct from graphene The STI phase contains 
unusual two-dimensional edge states, which may provide a route to spin-charge 
separation in higher dimensions [3, S] and may realize other novel magnetic 
and electronic properties [9]. It has been suggested that the interface of an STI 
and a conventional superconductor can support exotic quasiparticle states which 
could be used in fault-tolerant computing schemes |10[ |. It is currently believed 
that the Bii_a;Sba; alloys realize the only known topological insulator phase in 
the vicinity of a 3D Dirac point {s]. However, a particular challenge for the 
topological insulator Bii_a;Sba; system is that the insulating bulk gap is small and 
the material contains a significant amount of alloying disorder which is difficult 
to gate for the manipulation and control of charge carriers to realize a device. 
The topological insulator Bii_a^Sba; features five surface states of which only one 
carries the topological quantum number. Therefore, there is an extensive world- 
wide search for similar topological phases in stoichiometric materials with no 
alloying disorder, a larger gap and fewer yet odd numbered surface states. Here 
we present angle-resolved photoemission (ARPES) data and electronic structure 
calculation which seems to suggest the existence of a Z2 topological phase with a 
surface Berry's phase in the stoichiometric compound Bi2Se3. Our results taken 
collectively suggest that the undoped Bi2Se3 can serve as the parent matrix 
compound for the long-sought topological materials with only one surface state 
where in-plane carrier transport would be fully dominated by the Z2 topological 
effects. Moreover, the undoped compound should exhibit topological quantum 
effects at room temperature. 

Undoped Bi2Se3 is a semiconductor which belongs to the class of materials 312X3 
(X=S,Se,Te) with a rhombohedral crystal structure (space group Dg^ (i?3m)). The unit 
cell contains five atoms, with the quintuple layers ordered in the Se(l)-Bi-Se(2)-Bi-Se(l) 
fashion. These compounds have been studied extensively in connection to thermoelectric 
applications [12I. Experimental measurements suggest a band gap of approximately 0.35 
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where as the bulk band calculations estimate the gap to be in between 0.24 



to 0.32 eV pJl, Il7|]. Various measurements and calculations have reported that although 
the bulk of the material is a direct gap semiconductor [isl its electronic properties can vary 
significantly with different sample preparation conditions [3] , with a strong tendency to be 
n-type [18] due to impurity and vacancies. 

Single crystal of Bi2Se3 was grown by melting stoichiometric mixtures of high purity 
elemental Bi and Se in a 4 mm inner diameter quartz tube. The sample was cooled over a 
period of two days, from 850 to 650 °C, and then annealed at that temperature for a week. 
Single crystals were obtained and could be easily cleaved from the boule. High-resolution 
ARPES measurements were then performed using 17 to 45eV photons on beamline 5-4 at the 
Stanford Synchrotron Radiation Laboratory (SSRL) and beamline 12.0.1 of the Advanced 
Light Source at the Lawrence Berkeley National Laboratory. The energy and momentum 
resolutions were 15meV and 2% of the surface Brilloin Zone (BZ) respectively using a Scienta 
analyzer. The samples were cleaved in situ between lOK and 55K under pressures of less 
than 5 x 10~^^ torr, resulting in shiny fiat surfaces. For the deposition-efiFect measurements, 
the Fe atoms were deposited using an e-beam heated evaporator at a rate of approximately 
0.12A/min. We also carried our surface band calculations for the Bi2Se3(lll) surface for 
comparison with the experimental data. The calculations were performed with the LAPW 
method in slab geometry using the WIEN2K package {21]. GGA of Perdew, Burke, and 
Ernzerhof |22] was used to describe the exchange-correlation potential. SOC was included 
as a second variational step using scalar-relativistic eigenfunctions as basis after the initial 
calculation is converged to self-consistency. The surface was simulated by placing a slab of 
twelve quintuple layers in vacuum. A grid of 21 x 21 x 1 points were used in the calculations, 
equivalent to 48 k-points in the irreducible BZ and 300 k-points in the first BZ. 

Figure [I](b) shows the calculated surface electronic structure of the Bi2Se3(lll) along the 
K — T — M direction. The result with and without SOC are overlaid together for comparison. 
The bulk band projections are represented by the shaded areas in blue with SOC and green 
without. In the bulk, time-reversal symmetry demands E{k^ |) = E{—k^ [) while space 
inversion symmetry requires that E{k^]) = £"(— /c, |). Therefore, all the energy bands 
should be doubly-degenerate. However, space inversion symmetry is broken at the surface, 
so the degeneracy of the SS can be broken by spin-orbit interactions. Nevertheless, Kramer's 
Theorem dictates that spin degeneracy should be preserved at some high symmetry points 
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FIG. 1: Spin-orbit interaction induced surface Fermi surface: (a) A schematic of the bulk 3D BZ 
of Bi2Se3 and the 2D BZ of the projected (111) surface, (b) LDA band-structure of the 2D surface 
states along the K — T — M k-space cut. Bulk band projections are represented by the shaded 
areas. The band-structure results with spin-orbit coupling (SOC) is presented in blue and that 
without SOC is in green. No pure surface band is observed within the insulating gap without SOC 
(black lines). One pure gapless surface band crossing is observed when SOC is included (red 
lines), (c) The corresponding surface FS is a single ring centered at F. The band responsible for 
this ring is singly degenerate. The time-reversal-invariant momenta (TRIM) on the (111) surface 
BZ : the f and the three M points are marked by red dots. 

of the surface BZ that are time reversal invariant momenta (TRIM). In the 612863(1 11) 
surface BZ, these are given by F and three M points, located 60° away from each other and 
in the middle of two F points (Fig. [11(c)). The calculated result without SOC shows a clear 
gap between the valence and conduction bulk bands, with no pure surface bands crossing 
Ep. SOC drastically changes the band structure of the SS. One finds two singly degenerate 
surface bands emerging from the bulk projection which are paired together at the F point 
(one of the Kramers' points). The top surface band forms an electron pocket at E^?, giving a 
ring-like surface Fermi surface (FS). We emphasize that while doubly degenerate states may 
appear as two overlapping bands in the calculated band structure for some k- values, the two 
surface bands each appear only once, and represent one eigenstate at a given momentum 
and most importantly, the central Fermi surface is singly degenerate. 

The calculated valence band dispersion is in good correspondence while making a com- 
parison with the experimental data taken in both the T — M and T — K directions (Fig. [2j). 
The first and twelfth quintuple blocks are taken as the surfaces of the slab. The fraction of 
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FIG. 2: Experimental band dispersions are compared with theoretical calculations: The band 
second-derivative images and first-principles calculation results along (a) f — M and (h) T — K 
directions are presented. The color of the calculated bands represents the fraction of electronic 
charge residing in the surface layers. A rigid shift of E^^^ is included to match the lowest energy 
excitations in the ARPES data with calculations, a consequence of the system being somewhat 
electron-doped. The strongest signals are observed from the surface states (see text). 
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FIG. 3: kz dependence of low-lying excitations near f : High resolution surface band dispersion 
data near the f point along the (a) T-M and (b) T-K directions measured at 22eV photon energy 
are presented. The "V" shaped pure surface state (SS) band is nearly isotropic along both axes, 
forming a cone in the kx — ky space, (c) Results from the energy dependence study is presented 
along T-M going from F at 21 eV to 0.41 F — Z at 31 eV in leV intervals. Solid arcs indicate points 
in the bulk 3D BZ seen by the detector over a range of ±30°. (d) The energy dispersion data at 
three selected energies show no kz dependence for the outer SS band as well as the inner resonant 
state (RS) feature. At each energy the corresponding energy distribution curve is presented (top). 
The calculated bandwidth of the lowest 3D bulk conduction band (inset) is about 0.3 eV which is 
significantly larger than the instrumental resolution. 
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electron charge residing in the atomic spheres of these two surface blocks are presented by 
the color of the bands. Bright lines indicate bands located inside the bulk sample. These 
bands should overlap with a bulk energy band at a particular value, after projected onto 
the 2D surface BZ. A good agreement is seen by shifting the of the calculated result to 
20 meV above the bottom of the lowest conduction band - a consequence of doping the semi- 
conductor matrix with electrons. While this observation suggests that the bulk is slightly 
electron doped (n-type) as observed in many transport measurements, it is an advantage 
for the fact that we can have spectroscopic (ARPES) look at the large part of the complete 
band-structure of the surface states which was not possible in fully insulating BiSb A 
complete spectroscopic view is of significant advantage for we would like to study the con- 
nectivity of the surface and bulk states at all energies between the valence and conduction 
bands if possible to determine the unique and specific class of topological order of the parent 
materials j^. 

In the measured spectra, the strongest quasiparticle signals are typically observed near 
normal electron emission, at around 1.5-2.0 eV. This band corresponds to a pure surface 
band lying outside the bulk projection of states. Strong quasiparticle signals near T are of 
surface origin, based on their direct correspondence with the band calculation. Our energy 
dependent study of the bands suggest insignificant kz dependence of these state, confirming 
the 2D (surface) nature of the bands. Near E^?, one finds an additional 150 meV electron 
pocket-like feature inside the spin-split surface band. This feature is a hybrid resonant state 
formed by the superposition of the surface and the bulk. 

To investigate the resonant state (RS) feature centered at F, high momentum resolution 
data is presented in Fig. [3l The extracted band velocity of the outer pure SS in both 
T-K and F-M directions are approximately 5 x 10^ m/s, close to the calculated values. The 
SS crosses Ei. at 0.09 A"^ in the F-M direction and 0.11 A"^ in the T-K direction. The 
dependence of the bands at three selected energies and their respective energy distribution 
curves (EDCs) are presented in panel (d). By changing the energy of the incident photon, 
one can move to different k^ values in the 3D bulk BZ (Fig. [3](c)). The inner potential 
Vo used is approximately 10 eV, given by the muffin-tin zero of the calculation. The value 
is obtained by averaging the atomic potential over the interstitial region, using a muffin-tin 
radius of 2.5 bohr. Moving the photon energy at normal emission from 21eV, corresponding 
to F in the bulk 3D BZ, to 31eV, corresponding to 0.41 F — Z(or -0.59 T — Z after mapping 
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FIG. 4: The surface Fermi surface : (a) the surface FS is a small pocket around F. (b) High 
momentum resolution data around F show a single ring formed by the pure surface state band. 
In the middle of the ring is a filled-in disk-shaped spectral intensity reflecting resonant states (see 
text). The pure surface Fermi surface of BiSe is different from that of spin-orbit pair Fermi surface 
observed on metallic gold Au(lll). (c) The Au(lll) surface FS, which has two rings (each non- 
degenerate) surrounding the F TRIM. An electron circling the gold FS can only carry a Berry's 
phase of 0, characteristic of a trivial topological metal, Z2 ==+1. The single surface Fermi surface 
observed in Bi2Se3 reflects its non-trivial topological character Z2 =-1. 

into the first BZ), one finds no kz dispersion of the bands. A variation of the quasiparticle 
intensity is however observed with changing photon energy, due to matrix element effects. 
The calculated T — Z dispersion of the lowest conduction band is presented in the inset. The 
resulting bandwidth is approximately 300 meV. Therefore, if the features above 0.15eV are 
purely due to the bulk, one should observe a clear dispersion as moves away from the F 
point. The two closest bulk bands - the highest bulk valence band and the next conduction 
band, are at least 500 meV away with bandwidths of approximately 200 meV, so they are 
not to affect the quasiparticle signal near E^. This finding confirms that the outer "V" 
band as well as the inner electron pocket-like feature are due to the surface. 

The surface FS is then presented in Figure [H Panel (a) presents the ARPES data 
over the entire 2D (111) surface BZ. No E^? crossings are observed other than the features 
centered at F. The three TRIMs located at M are not enclosed by any Fermi arcs. This 
observation is in contrast to the surface FS of Bii-^^Sb^; [3], which contains a electron-hole 
pocket features around each of the M points. The region around F is zoomed in with 
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finer momentum resolution data (panel (b)). One finds a ring-like feature formed by the 
outer "V" SS surrounding the filled-in RS disk. This ring is single degenerate thus spin- 
polarized as expected from the correspondence bewteen the LDA band calculation. An 
electron encircling the surface FS that encloses a TRIM collects a Berry's phase of tt mod 
27r, picking up a value of tt from each of the Fermi arcs that enclose the TRIM. At the 
first sight the surface FS in Bi2Se3 may look like the well-known surface FS in bulk metallic 
strong spin-orbit coupled gold (Au). However, the SS in gold(lll) is split into two parabolic 
bands each singly degenerate, which are shifted in momentum-space from each other and 
both enclose the F-point {24, 25]. The resulting surface FS topology contains two concentric 
rings around F, instead of one as observed in our data for Bi2Se3. Therefore, the Fermi arcs 
in gold encloses the TRIM an even number of times, with two crossings between a pair 
F and M points. Because each of the rings is singly degenerate (necessarily spin polarized), 
an electron encircling the FS can only carry a trivial Berry's phase of mod 27r, by picking 
up one TT from each of the rings. This makes gold surface states topologically trivial. 

To determine the topological nature of Bi2Se3, one needs to determine if the RS con- 
tributes an odd or even number of crossings between two TRIMs. In order to determine 
the nature of the RS states and its connection to the conduction band structure we need 
to doped the surface of Bi2Se3 with atoms that donates electrons on the surface. This can 
be done by depositing alkali atoms on the cleaved surface under UHV conditions. For this 
particular experiment, we have chosen to deposit iron on the surface. Figure [5](a) presents 
ARPES data at 29eV along T — M after Fe atoms of less than 1% of a monolayer is deposited 
onto the sample. The system becomes strongly electron doped, with the bottom of the "V" 
band shifted by more than 200meV. In addition to a shift in the chemical potential, the 
outer "V" band also becomes less dispersive, with a decreases in the band velocity near Ep. 
More interestingly, the inner RS appears to be the bottom of a pair of spin-split parabolic 
bands shifted from each other in k-space. These bands are reminiscent of those observed 
in the SS of Au(lll) [24, M and Sb(lll) [27] at F. The result shows that the RS feature 
crosses Ep^ twice, bringing the total number of crossings between a T-M pair to three. 

Since iron is magnetic its use in doping the surface also serves a separate but important 
purpose which has to do with the possiblity of observing the system's response to the break- 
ing of time reversal invariance on the surface. With magnetic irons on the surface Kramer's 
degeneracy is broken at the TRIMs and a gap is expected to open up. A small gap or 
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FIG. 5: The iron deposited surface: band structure in a broken time-reversal setting and the 
determination of contribution of the RS feature to the Berry's phase counting: (a) a pair of spin- 
spht bands are extended from the RS feature after the chemical potential is shifted with doping. A 
small gap opens at F, due to the loss of time reversal symmetry. The corresponding EDC is shown 
in (b). Similar to (c) n-Bi2Se3, (d) the band structure of the strong topological metal Sb(lll) 
contains RS at some TRIMs. The pair of spin-split RS bands in Bi2Se3 resembles the SS in Sb. In 
both cases (e) a "partner-switching" behavior is observed, a critical characteristics in generating 
an odd number of surface state Fermi crossings. 

strong spectral weight suppression is indeed observed at f (Fig. [5](a)). The outer "V" band 
becomes more parabolic and is detached from the lower "A" band. The signal intensity at F 
becomes very weak, whereas in the undoped case, the signal is strongest at the Dirac point. 
For this reason, an exact measurement of the gap size is difficult. Spectral suppression or 
a weak gap opening at the Kramers' point suggest that magnetic field and disorder can 
cut the k-space surface band thread that connects the bulk valence and conduction bands 
driving the topological insulator into a trivial band insulator as expected in Z2 theory of 
these materials |26l | 

Therefore, the FS of the Bi2Se3(lll) SS only encloses the F point an odd number of 
times, giving it a z/q = 1 and Z2=-l band topology similar to Bii_a;Sba;. However, the 
surface transport of pure undoped Bi2Se3 would then be dominated by topological effects 
since there is only one surface band that is also topological. Our calculation suggests that 
the topological character should be preserved in the undoped compound, which is insulating 
with metallic boundary states. One might be able to obtain the pure undoped compound 
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by hole-doping the n-type Bi2Se3 thereby shifting the chemical potential downwards. Such 
a chemical alteration would also remove the resonant state discussed previously, leaving a 
single ring in the surface FS which carries a tt Berry's phase. Treated along these lines our 
data in Fig. -3 and 5 suggest that the undoped material should have a band gap of about 0.3 
eV which is sufficient to keep it insulating at room temperature. 

In conclusion, we have calculated the band structure of Bi2Se3(lll) surface and found that 
spin-orbit coupling induces a single non-degenerate band crossing at E;?. Our experimental 
data agrees with our surface band calculations as far as the band topology is concerned which 
demonstrates that the stoichiometric compound belongs to the Z2—-I topological class. We 
have also shown that since that Fermi surface contains only one surface state, the transport 
properties of undoped 612863 would be dominated by topological effects and since the band 
gap is large (Fig.-3 and 5) the material can be considered as a room temperature topological 
insulator without any alloying disorder. 
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